Penta-Hepta Defect Chaos in a Model for Rotating Hexagonal Convection 
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In a model for rotating non-Boussinesq convection with mean flow we identify a regime of spatio- 
temporal chaos that is based on a hexagonal planform and is sustained by the induced nucleation 
of dislocations by penta-hepta defects. The probability distribution function for the number of 
defects deviates substantially from the usually observed Poisson-type distribution. It implies strong 
correlations between the defects in the form of density-dependent creation and annihilation rates 
of defects. We extract these rates from the distribution function and also directly from the defect 
dynamics. 
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Spatio-temporal chaos is at the focus of experimental 

SlJ J, i i I, I i and of theoretical 0, 0, 
]l4 , |l5| research in high-dimensional dynamical sys- 
tems. Most of the extensive studies have been devoted 
to variants of thermally or electrically driven convection 
in thin liquid layers ^, ^ ||, |ll|, 14 . Detailed exper- 



imental studies have also been performed on vertically 
vibrated layers of fluids ^ and on chemical systems [||. 
Theoretically, various regimes of spatio-temporal chaos 
of the complex Ginzburg-Landau equation have been in- 
vestigated g |l§. 

A striking feature of most spatio-temporally chaotic 
states are defects in the pattern. They can be line defects 
like domain walls, point defects like dislocations, disclina- 
tions, and spirals, or composite defects like penta-hepta 
defects. In particular dislocations have attracted great 
attention since they are easy to identify. Investigators 
have utilized their statistical, geometrical and dynami- 
cal aspects to quantify the chaotic states in which they 
arise. For example, the number of dislocations (spirals) 
in the wave patterns governed by the complex Ginzburg- 
Landau equation has been found to obey Poisson-type 
statistics This suggests the interpretation that in 

this system dislocations are created randomly in pairs 
with a fixed probability, after which they diffuse through- 
out the system without any mutual interaction until they 
annihilate each other in collisions [|l^. The correspond- 
ing behavior and associated distribution function have 
also been found experimentally in electrically driven con- 
vection in nematic liquid crystals |l|] and in thermally 
driven convection in an inclinded layer and theoret- 
ically in simulations of coupled Ginzburg-Landau equa- 
tions for parametrically excited standing waves |l5| . 

Geometric aspects of dislocations have been investi- 
gated in experiments on binary- mixture convection where 
the possibility to reconstruct the patterns from the dis- 
locations has been explored [^. In another study the ge- 
ometry and connectivity of the dislocations' world lines 
in space-time has been considered . Through the cre- 
ation and annihilation events the world lines form loops 
in space-time. In studies of a type of defect-unbinding 
transition it has been found that the degree of order of 



the defected pattern is related to the statistics of the size 
of the loops. 

The dynamical relevance of dislocations has been sug- 
gested in direct simulations of the Navier- Stokes equa- 
tions of spiral-defect chaos in Raleigh-Benard convection. 
It was found that the chaotic state is by far most sensitive 
to p erturbations during the creation of dislocation pairs 
[|l4[ . The best evidence for the significance of defects as 
dynamical objects has been provided in simulations of 
the complex Ginzburg-Landau equation where the con- 
tribution of the defects to the Lyapunov dimension of the 
chaotic attractor has been extracted 

Most of the detailed analyses of spatio-temporal chaos 
and of its defects have been performed in disordered pat- 
terns that are based on stripes (or rolls). Much less work 
has been done on spatio-temporal chaos related to other 
planforms like rectangles Jo] (and, related to it, vector 
waves 0) or hexagons [|[, and the role of the corre- 
sponding defects has been barely touched upon. 

In this Letter we describe a spatio-temporally chaotic 
state that is based on a hexagonal pattern. Its disorder 
is closely tied in with the appearance of penta-hepta de- 
fects (PHDs), each of which consists of two dislocations 
in two of the three modes making up the hexagon pat- 
tern. In contrast to most other systems discussed above 
it is not only the instability of the background pattern 
that drives the chaotic state, but also the instability of 
the PHDs themselves. Thus, in the presence of PHDs 
new dislocations are created through induced nucleation. 
As a consequence the probability distribution function 
for the number of defects is considerably broader than 
the Poisson-type distributions reported in previous stud- 
ies ||, |lO|. We obtain this persistent, chaotic state in a 
Swift-Hohenberg-type model for rotating non-Boussinesq 
convection at low Prandtl numbers. While induced nucle- 
ation itself has been reported previously |l^ , without 
rotation it did not sustain persistent chaotic dynamics 

Motivated by the strong effect of mean flows and ro- 
tation on convection roll patterns ||, || we have previ- 
ously studied their effect on the stability of hexagon pat- 
terns and their PHDs within the framework of Ginzburg- 
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FIG. 1: a) Snapshot of penta-hepta defect chaos for a = 0.4, 
7 = 2, ;5 = -2.6, R = 0.17, and L = 233. b) Corresponding 
space-time diagram of the radially integrated Fourier spec- 
trum. A movie of the temporal evolution of the pattern and 
of its defects can be found on the EPAPS server. 

Landau equations [|l9|, Since the Ginzburg-Landau 
equations break the isotropy of the system they are not 
suited for investigations of spatio-temporal chaos. In this 
paper we therefore investigate a minimal extension of the 
Swift-Hohenberg model, 

dt^ = R^- (V' + if^ - + a (V^)' + (1) 

The quadratic terms proportional to a and 7 break 
the up-down symmetry ip ^ —ip and model the non- 
Boussinesq effects. The chiral symmetry is broken by 
the terms involving 7 and J; thus, to leading order these 
coefficients are linear in the rotation rate. The mean- 
flow velocity and its stream function are given by U and 
^, respectively, and 13 increases with decreasing Prandtl 
number. We simulate (|^) numerically using a parallel 
pseudospectral code with periodic boundary conditions. 
A typical snapshot in the chaotic regime (Figj^a) shows 
domains of hexagons of distinct orientations separated by 
domain walls in which many PHDs are aggregated. Due 
to the broken chiral symmetry most of the domains pre- 
cess slowly counterclockwise. The corresponding space- 
time diagram for the temporal evolution of the radially 
integrated Fourier spectrum is presented in Fig.p^. 

To identify the dislocations and PHDs we make use of 
the fact that despite the disorder of the pattern its spec- 
trum exhibits six peaks that are clearly separated most 
of the time and that are rotated by 120^ with respect to 
each other (cf. Fig.p^). We demodulate the pattern using 
three carrier wavevectors that slowly precess along with 
the spectrum, ip = Y^^=i exp{iq^j{t) • r) + h.o.t. + c.c. 
Figs.||a-c show the temporal evolution of a smaller section 
of the pattern with the dislocations in the three modes 
marked by triangles, squares, and circles, respectively. 
Open (closed) circles denote a positive (negative) topo- 
logical charge of the dislocations. 



FIG. 2: Induced nucleation of dislocations. Enlargements 
corresponding to the box in Figjl| at times t = 747, t = 759 
and t = 760. Dislocations in Ai, i = 1..3, marked by squares, 
circles, and triangles, respectively. Open (filled) symbols in- 
dicate positive (negative) topological charge. White (grey) 
lines mark heptagons (pentagons) making up the PHDs. 
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FIG. 3: Probability distribution function for the number of 
dislocation pairs in the pattern, the parameters are as in Figjl| 
for L 233. For the inset L 114, a 0.4, 7 = 3, /3 = -5 
and R = 0.09. Solid line is fit to (§). 



In various experimental and theoretical investigations 
of stripe-based disordered patterns the probability distri- 
bution function for the number of defects has been used 
to obtain a first characterization of the defect evolution 
^ H, |l^. Except for the ordered chaotic state in 
the probability distribution function for the number of 
defects were found to be close to a Poisson-type distribu- 
tion, indicating that the dynamics are consistent with the 
simple diffusive model described above with very weak 
correlations between the defects [|l^. In particular, the 
creation rates depend only little on the defect density 
[||. However, this is not the case for the defect chaos in 
hexagons. Fig.|| gives the distribution function for the 
number of dislocations in the penta-hepta defect chaos 
for two system sizes, L = 233 and L = 114 (inset), and 
two sets of parameter values. The symbols give the rel- 
ative frequency to find n dislocation pairs in one of the 
three modes, whereas the dashed line gives the best fit to 
the squared Poisson distribution (with the same mean) 
corresponding to the uncorrelated dislocation dynamics 
[|lO|. Clearly, in the penta-hepta defect chaos the defect 
dynamics are far from uncorrelated. 

A more detailed analysis of the defect dynamics re- 
veals a strong tendency for dislocations to be created in 
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FIG. 4: Persistence limit for penta-hepta defect chaos as a 
function of p and 7. Other parameters as in inset of Figjs] 
(L = 114). 



the vicinity of already existing PHDs. This is illustrated 
in Fig.||. Due to the gradient terms involving a and 7, 
which lead to nonlinear gradient terms in the Ginzburg- 
Landau equations the dislocations making up the 
PHDs are spatially separated [|l^, |9| (cf. Fig.^). In 
addition, a PHD in modes Ai and say, leads to a 
perturbation in mode As. For sufficiently large a and 7 
the perturbation evolves into a dislocation pair in mode 
As (in Fig.||b splitting of the cell between the 'square' 
and the 'circle' dislocation). The newly created disloca- 
tions then recombine with the oppositely charged dislo- 
cations in the oris; inal PHD to form two PHDs (Fig.||c), 
which then typically move apart from each other. Such 
induced defect nucleation has been found previously in 
coupled Ginzburg-Landau equations |l^, ^ and in a 
Swift-Hohenberg-type model without rotation or mean 
flow Q. However, in contrast to the case discussed in 
[|l8[ , in the presence of rotation the nucleation is sufficient 
to sustain a precessing chaotic state. As shown in Figj^, 
for small Prandtl numbers {P < 0) mean flow enhances 
the persistence of the chaotic state. 

To establish a quantitative connection between the in- 
duced defect nucleation and the defect distribution func- 
tion we consider an extension of the simple kinetic model 
for the dynamics of the defects presented in [|l^. Since 
there are three different modes Aj and because the to- 
tal topological charge of a PHD has to vanish [|2j, the 
statistics of the defect dynamics are described by a com- 
bined distribution function (^12' ^12^ ^235 ^235 ^31 5 ^31) 
for the six different kinds of PHDs. Here denotes, 
for instance, the number of PHDs involving a disloca- 
tion with positive charge in mode Ai and a dislocation 
with negative charge in mode A2. In principle, there 
are also dislocations that are not bound in a PHD. In 
this kinetic model we assume that their dynamics are 
fast enough to follow quickly the number of PHDs. The 
change in Vq during a time interval At can be expressed 
in terms of creation and annihilation rates, which de- 
pend on the various defect densities. The numerical sim- 
ulations show that the densities are strongly correlated 
at equal times [pll, which implies that Vq is strongly 



peaked when its six arguments are equal. Integrating 
out the dependence of Vq on five of its arguments, one 
therefore obtains a closed approximation for the change 

in V{n^2 = ^) = / 7^6 <^^r2^^23"-^^3i during a time in- 
terval At involving the creation and annihilation rates 

^ n 1 

V{t + At, n) = V{t, n) + At {r+_^V{t, n - 1)+ 

r-+iP(t,n + i)-(r- + r+)p(t,n)}. (2) 

In steady state the distribution function satisfies detailed 
balance, V{n + 1) F"^^ = V{n)T^. Assuming a fixed 
rate for the induced nucleation, the rate for the process 
shown in Fig.|| depends linearly on the density of the 
'square-circle' PHDs. It creates one 'triangle-circle' (and 
one 'triangle-square') PHD and annihilates the original 
'square-circle' PHD. This suggests a linear contribution 
to the dependence of the annihilation and creation rates 
on the defect density. The reverse process originates from 
two PHD's and therefore contributes quadratic terms. 
Including also the spontaneous creation of dislocations, 
which then form PHDs, we make the ansatz 



Co + cin + C2n . 



(3) 



Since the probability distribution depends only on the 
relative rates we choose the overall time scale to normal- 
ize the coefficient a2 of the quadratic annihilation rate to 
unity. The steady-state solution to (0,0) is then given by 



Vin) = V{0) n 



j=0 



Co + CiJ + C2j" 

ai{j + 1) + {j + If 



(4) 



with 7^(0) determined by the normalization condition. A 
fit of the numerical simulation results to is shown as 
solid line in Fig.|| and its inset. For both system sizes 
the fits are very good. For L = 114 we obtain cq = 20.7, 
ci = 20.7, C2 = 0.12, and ai = 8.6 (a2 is scaled to unity), 
confirming the strong dependence of the creation rate on 
the number of defects. 

By tracking each dislocation from its creation to its 
annihilation we can also determine the creation and an- 
nihilation rates directly from the numerical simulations. 
Fig.^ shows these rates for a dislocation in a given mode 
as a function of the number of dislocation pairs in the 
same mode for a system of size L = 114 (same parame- 
ters as in inset of Fig.||). In principle, the rates should 
be given as functions of the number of PHDs involving 
the other modes. However, due to the finite distance 
between the dislocations within a PHD the grouping of 
dislocations into PHDs is not always unique. Because the 
numbers of dislocations in the three modes are strongly 
correlated, taking the number of dislocations in the same 
mode provides a good approximation. The large scatter 
in the data for larger defect numbers is due to the lack of 
statistics for events of that kind (cf. inset of Fig]|). Sim- 
ilarly, there are only few events with few defects. Clearly, 
in the intermediate range of n not only the annihilation 
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FIG. 5: Creation (squares) and annihilation (circles) rates of 
dislocations as a function of the number of dislocation pairs 
in the same mode. Parameters as in inset of Fig|3| (L = 114). 
Inset: wavenumber distribution function with stability limit 
for hexagons according to weakly nonlinear theory (dashed) . 

rate but also the creation rate depends strongly on the 
defect number. 

To connect the directly measured rates with the dis- 
tribution function (^) the solid curves in Fig.|| give the 
creation and annihilation rates as determined from fitting 
the distribution function for the defect number using the 
form (H). For this comparison the overall time scale is 
adjusted to fit the time scale of the simulations. Clearly, 
the rates inferred from Fig.^ agree quite well with the di- 
rectly measured ones over the statistically reliable range 
of n and confirm the interpretation of the deviation of 
the distribution function from the squared Poisson dis- 
tribution. 

The creation rate for dislocations does not vanish for 
n = 0, i.e. when no PHDs are present. This indicates the 
spontaneous creation of dislocation pairs directly from an 
instability of the hexagonal pattern, although the mean 
wavenumber of the background pattern is clearly inside 



the stability balloon. However, the distribution function 
for the local wavenumber (inset of FigjH]) shows that there 
is a noticeable tail of the distribution function that ex- 
tends beyond the low-q stability limit, as determined by 
a weakly nonlinear analysis of (|^). This suggests that 
some dislocation pairs are created through a side-band 
instability of the periodic pattern. 

In conclusion, in a model for rotating non-Boussinesq 
convection we have identified a spatio-temporally chaotic 
state that is dominated by the dynamics of penta-hepta 
defects of the underlying hexagon pattern. In contrast 
to previously analyzed chaotic states, which are stripe- 
based, the defect statistics of this penta-hepta chaos indi- 
cate strong correlations between the defects. We identify 
the origin of the correlations as the induced nucleation 
of dislocations due to the presence of penta-hepta de- 
fects. From the defect statistics we extract the depen- 
dence of the creation and annihilation rates of defects on 
the defect density and find good agreement with the rates 
measured directly by following individual defects in the 
simulations. In ongoing direct simulations of the Navier- 
Stokes equations for rotating non-Boussinesq convection 
we have identified regimes exhibiting induced nucleation 
of dislocations [^. In the simulations performed so far 
the induced nucleation either occurs only as a transient 
and eventually leads to ordered hexagon patterns or it 
leads to persistent chaotic dynamics that are somewhat 
more complex than the dynamics found here for the ex- 
tended Swift-Hohenberg model. 
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